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ABSTRACT 

We report on Australia Telescope observations of the massive star forming region 
G305.2+0.2 at 1.2 cm. We detected emission in five molecules towards G305A, confirm- 
ing its hot core nature. We determined a rotational temperature of 26 K for methanol. 
A non-LTE excitation calculation suggests a kinematic temperature of order 200 K. A 
time dependent chemical model is also used to model the gas phase chemistry of the 
hot core associated with G305A. A comparison with the observations suggest an age 
of between 2 x 10"* and 1.5 x 10^ years. We also report on a feature to the SE of G305A 
which may show weak Class I methanol maser emission in the line at 24.933 GHz. The 
more evolved source G305B does not show emission in any of the line tracers, but 
strong Class I methanol maser emission at 24.933 GHz is found 3" to the east. Radio 
continuum emission at 18.496 GHz is detected towards two Hii regions. The implica- 
tions of the non-detection of radio continuum emission toward G305A and G305B are 
also discussed. 

Key words: masers ~ stars: formation - infrared: ISM ~ ISM: molecules 



X 



1 INTRODUCTION 

The study of the early stages of massive star formation 
(MSF) in our Galaxy has, until recently, received relatively 
little attention compared to low mass star formation. How- 
ever, with the advent of telescopes, particularly in the mil- 
limetre (mm) and sub-mm part of the electromagnetic spec- 
trum, it has become possible to study the earliest stages of 
MSF. This is the strongest part of their spectral energy dis- 
tribution (SED) measureable from the ground. It also con- 
tains a myriad of complex molecular transitions, allowing 
investigation of the physics and chemistry surrounding these 
young sources. 

G305.2-I-0.2 is a site of MSF in the southern Galactic 
plane. Cl ass II methanol ma sers were reported in two po- 
sitions bv iNorris et~al] (|l993l ): G305. 21-^0.21 (G305A) and 
G305. 20-1-0.21 (G3 05B), separated by approximately 22". 
IWalsh et al] (|l999l ) observed this region in the near-infrared 
and found that whilst G305A is not associated with any in- 
frared source, G305B is coincident with a bright and very 
reddened inf rared source. Mo re sensitive near-infrared ob- 
servations bv lDe Buized l|2003l ) indicate that there is indeed 
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a weak infrared source coincident with G305A, however, it is 
not clear if it is associated with the maser site or is instead an 
unrelated foregroun d star. Obser vations in the mid-infrared 
(10.5 and 20 fim) bv lWalsh erall ((20011 ') confirm the infrared 
source associated with G305B has a steeply rising spectral 
energy distribution towards longer wavelengths. However, 
no mid-infrared source was found associated with G305A. 
The Galactic Legacy Infrared Mid-Plane Survey Extraordi- 
naire (GLIMPSE) is a legacy science program of the Spitzer 
Space Telescope. It confirms the detection of a weak infrared 
source at the position of G305A at 3.6 and 4.5/^m. However, 
the region shows much extended emission (see Figures la 
and 3b), which makes photometry unreliable. 

et al. I (|2005l ) have detected a strong 1.2 mm contin- 
uum souTce_whflse_ge^c_coincides with G305A (see figure 
1 of I Walsh fc BurtonI l|2006f) '). The 1.2 mm map also shows 
some extended emission in the direction of G305B, but the 
extended emission was poorly resolved. Because G305B ap- 
pears to dominate the luminosity at shorter wavelengths (up 
to about 20 fim) and G305A appears to dominate at longer 
wavelengths (1.2 mm), there will be a transition point at far- 
infrared wavelengths where the two sources will have equal 
brightnesses. Thus, it is impossible to tell which source is the 
most luminous, based on infrared data alone. Neither maser 



2 A. J. Walsh et al. 



site is coincident with any detected radio continuum source 
l|Phillips et al ]|l998al ). However, one mig ht be expected from 
G305B, since extrapolation of its SED ijWalsh et al.ll200lh 
indicates it is powered by a star bright enough to produce an 
observable ultracompact (UC) Hii region. Furthermore, it is 
not known why the brightest maser site, G305A is not asso- 
ciated with a bright infrared source, nor a radio continuum 
source and yet appears to be a strong source of continuum 

emi ssion at 1.2 mm. 

IWalsh fc BurtonI (|2006l ) found G305A to be the local 
peak of emission for 1-0 transitions of ^^CO, HCO^ and 
N2H+, as well as transitions of CH3CN and CH3OH. NH3 
(1,1), (2,2) and (4,4) emission has been detected coinci- 
dent with G305A (Steven Longmore, private communica- 
tion 2007). This suggests G305A is a hot core and may be 
the most massive component in the region. The line and 
dust continuum emission appear to be slightly elongated in 
the direction of G305B. Therefore it appears G305B may be 
associated with a secondary hot core, wh ich was spatially 
unreso lved with the previous observations. I Walsh fc BurtonI 
also identified the possible prestellar massive core 
G305SW, which was only detected in N2H"'^ and shows nar- 
row lines, compared to the emission elsewhere in the field 
of view. In this paper we attempt to better understand the 
nature of these sources by investigating both molecular line 
and continuum emission at 1.2 cm. 



2 OBSERVATIONS 

Observations towards G305.2-f0.2 were carried out using 
two configurations of the ATCA in 2005. The first observa- 
tions took place on 2005, May, 3-4th and utilised the 750A 
configuration. The second observations took place on 2005, 
September, 4th and utilised the H214 configuration. Details 
of each observing session are given below. 

2.1 2005, May observations 

Using the 750A configuration, we observed G305 with base- 
lines ranging from 77 to 735 m, with five antennas, and base- 
lines up to 3750 m when including antenna 6. However, data 
from antenna 6 were only used when observing masers. The 
primary beam for these observations was between 2.0' and 
2.6', and the synthesised beam was between 3.2" and 4.9", 
when only 5 antennas were used. Observations were made to- 
wards a single position, with pointing centre 13** 11"* 10.6", 
right ascension and -62° 34' 38" declination (J2000). We 
spent a total of 12 hours observing G305 in this configura- 
tion, which included time for primary calibration on 1934- 
638, bandpass calibration on 1253-055 and phase calibration 
on 1414-59. We observed in a frequency-switching mode, 
with main (spectral line) frequencies being 18199 MHz and 
24932 MHz, which include emission lines of cyanoacetylene 
(HC3N (2-1)) and the K=2 ladder of methanol (CH3OH), 
respectively. Assumed rest frequencies for these lines are 
given in Table [1] Each frequency was observed with 256 
channels over a bandwidth of 16 MHz, resulting in a velocity 
resolution of l.Okms"^ and 0.75 kms~^ per channel at the 
HC3N and CH3OH line frequencies, respectively. In addi- 
tion to the spectral lines, continuum emission was observed 
at 18496 MHz with a bandwidth of 128 MHz. 



2.2 2005, September observations 

Using the H214C hybrid configuration, we observed G305 
with baselines ranging from 82 m to 240 m with five an- 
tennas. The primary beam for these observations was be- 
tween 2.0' and 2.6', and the synthesised beam was be- 
tween 8.6" and 11.8". Observations were made with a sin- 
gle point, with pointing centre 13'' 11"* 6.25", right as- 
cension and -62° 34' 53" declination (J2000). The reason 
for the different pointing centre to the 2005, May obser- 
vations was to cover th e newly discovered source G305SW 
(j Walsh fc BurtonI I2OO6I ). We spent a total of 11 hours ob- 
serving G305 in this configuration, which included time for 
primary calibration on 1934-638, bandpass calibration on 
0537-441 and phase calibration on 1414-59. We observed in 
a frequency-switching mode, with five pairs of spectral line 
frequencies, the details of which are given in Table [1] 

Observations made in this session have poorer spatial 
and spectral resolution than the 2005, May session. The 
main implication of this is that the spectral resolution is in- 
sufficient to identify most systemic motions, although we are 
able to identify detections and discriminate between emis- 
sion from G305A and G305B. 



3 RESULTS 

Approximately half the lines that were searched were de- 
tected, and details of these lines are given in 

TableO The ta- 
ble also includes details of those lines that were not detected, 
together with their upper limits. Since the 2005 September 
observing run was conducted using a relatively poor spa- 
tial and spectral resolution, we do not know line widths or 
precise locations for most of the lines detected during this 
time. In the last column of Table [2] we list the position of 
the emission and nearest association if it is not found coin- 
cident with G305A or G305B. The observations of methanol 
and cyanoacetylene conducted in 2005 May have both higher 
spatial and spectral resolution, and so we report more de- 
tailed results on each below. 



3.1 CH3OH (methanol) 

The bandpass covered three lines of methanol during these 
observations, as listed in Table[T] We find strong, narrow line 
emission from the J=4 line close to the position of G305B, 
but not coincident (Figure [1]) . The emission appears to be 
offset from G305B by approximately 3" to the east (13 11 
10.89 -62 34 38.5 J2000). In addition to this, weaker emis- 
sion is seen in all three methanol lines at the position of 
G305A (Figure [1]). The methanol emission around G305A 
also shows extended emission to the south-south-east, ex- 
tending by approximately 5". This extension, which we refer 
to as G305A(SE), is strongest in the J=4 methanol line, but 
the other two lines are also detected here. We note also that 
the J=4 methanol line appears to be very narrow at this 
position. So the relative methanol line intensities and line 
width are reminiscent of a weaker version of the emission 
seen close to G305B. 
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Table 1. Properties of observed lines. The notation for each line in the second column either conforms to the structure of (Juppor - 

Jlower) or (J-K, + K - J-K,+k)- 



Date 


Observed 


riest 
Fr e G ncy 
(MHz) 


Primary 
Beam 
(') 


Synthesized 
Beam 
(") 


Bandwidth 
I iviinz ) 


velocity 
Resolution 

(kms"-'-) 


2005 MAY 3-4 


HC3N (2-1) 


18196.226 


2.6 


4.9 


16 


1.0 


2005 MAY 3-4 


CH3OH (32,1 - 3i,2) E 


24928.700 


2.0 


3.2 


16 


0.75 




CH3OH (42,2 - 4i.3) E 


24933.470 












CH3OH (22,0 - 2i,i) E 


24934.380 










2005 SEP 4 


NH3 (1,1) 


23694.496 


2.0 


9.0 


32 


3.2 




NH3 (2,2) 


23722.633 


2.0 


9.0 


32 


3.2 


2005 SEP 4 


H2O (61,6 - 52,3) 


22235.080 


2.1 


9.6 


32 


3.4 




CH3OCHO (2i,2 - li,i) E 


22827.741 


2.1 


9.4 


32 


3.3 




CH3OCHO (2i,2 - li,i) A 


22828.134 


2.1 


9.4 


32 


3.3 




NH3 (3,2) 


22834.185 


2.1 


9.4 


32 


3.3 


2005 SEP 4 


CH3OH (82,6 - 81,7) E 


25294.417 


1.9 


8.6 


64 


11.9 




DC3N (3-2) 


25329.437 


1.9 


8.5 


64 


12.0 


2005 SEP 4 


OCS (2-1) 


24325.927 


2.0 


8.8 


64 


12.3 


2005 SEP 4 


CH3CN (1-0) 


18399.892 


2.6 


11.6 


32 


4.1 


2005 SEP 4 


DC3N (2-1) 


16886.310 


2.8 


12.7 


64 


17.8 




H91/3 (RRL) 


16894.200 


2.8 


12.7 


64 


17.8 




HC7N (15-14) 


16919.979 


2.8 


12.7 


64 


17.7 




HC5N (7-6) 


18638.616 


2.6 


11.5 


32 


4.0 



3.2 HC3N (cyanoacetylene) 

We only detect cyanoacetylene coincident with the po- 
sition of G305A (Figure [2}. The emission appears 
to be unresolved spatially, but is resolved spectrally 
with a line width of 5.1km s~^. Cyanoacetylene can be 
formed by ne utral-neutral reactions in interstellar space 
l|Fukuzawa fc O samura 1997^ ), and so is likely a precursor 
molecule th a t is p resent at the onset of star formation. 
iMillar et al.l |l993) indicate that cyanoacetylene can also 
be formed in hot cores. Their models also suggest that the 
cyanoacetylene abundance is greatly reduced after about 10^ 
years. Implications of the time-dependant chemistry are dis- 
cussed further in fjSl 



3.3 Continuum emission 

Our 2005 May data included a frequency (18496 MHz) de- 
voted to continuum emission, with a bandwidth of 128 MHz. 
The continuum map is shown in Figure [S] Two promi- 
nent continuum fea tures appear in the m ap: the compact 
Hii region G305HH (jWalsh fc Burtonll2006l ). and a similarly 
sized compact Hii region about 30" to the SE of G305A, 
which we call G305HII(SE). We assume that each Hii re- 
gion has its own powering source that is unrelated to ei- 
ther G305A or G305B. This is because there is a large 
separation between the Hii regions and G305A or G305B. 
Also, the powering source for G305HH appears as a near- 
infrared source at the centre of the radio continuum con- 



tours J Walsh fc Burtonll2006l ). The integrated intensity of 
the emission from G305HH is 174 mjy and the integrated 
intensity of the emission from G305HH(SE) is 57mJy. Line- 
free channels from the 2005 September data can be used 
to reconstruct the fluxes of the Hii regions between 16 and 
25 GHz. A plot of the distribution of fluxes with frequency is 
shown in Figure ID Overall, it appears that the flux remains 
constant over the frequency range, which is consistent with 
optically thin free-free emission. There are apparent changes 
in flux in both sources which are matched in frequency. This 
is most likely due to errors in the absolute flux calibration, 
which precludes us from ass essing small changes i n the slope 
of the continuum emission. IPhil lips et all (|l998al ) have pre- 
viously measured the continuum flux of G305HH at 8.6 GHz 
to be 195 mJy, which is also consistent with optically thin 
free-free emission. 

We can make an estimate of the type of stars that 
give rise to these Hii regions, using equation (4) of 
iKurtz. Churchwell fc Wo"odl (|l994l ). Here, we assume the dis- 
tance to G305 is 3.9 kpc, the electron temperature (Te) 
is IQQQQK and a(i/, Te) — 1. If we also assume that no 
photons are absorbed by dust within the Hii regions, then 
we find the rate of emission of lyman continuum pho- 
tons to be 2.7 X lO^'^s"^ and 8.8 x lO^^s'^ for G305HII 
and G 305HII(SE), respectively. Using the tables of iPanagial 
(|l973t) . these rates are equivalent to slightly earlier than a 
BO spectral type for G305HII and slightly later than a BO 
spectral type for G305HII(SE). If, however, 90% of the ly- 
man photons are absorbed by dust within the Hii regions, 



4 A. J. Walsh et al. 



CHgOH (4gg - 4i3) E Integrated Intensity 



CHgOH spectrum at G305B 



-63''34'30" - 




13''11™14= 13^ 12^ 11^ 

Right Ascension (J2000) 




3,, - 3,„ 



Velocity (km/s) 

CligO;] =-pecl.rurn al G305A(SK) 




Velocity (km/s) 



Figure 1. (a) CH3OH integrated intensity map of the J=4 line. Plus symbols show the positions of the 6.7 GHz CH3OH maser sites 

G305A and G305B. Contour levels arc 0.2, 0.4 and 0.6Jy/bcam.km/s. (b) CH3OH spectrum towards the bright line 3'' to the cast of 
G305B, at 13 11 10.89 -62 34 38.5 (J2000). (c) CH3OH spectrum towards G305A. (d) CH3OH spectrum towards G305A(SE), at 13 11 
14.00 -62 34 46.5 (J2000). 



HC3N (3-1) Integrated Intensity 




Right Ascension (J2000) 



HC3N (2-1) spectrum ot G305A 




-50 
Velocity (km/s) 



Figure 2. (a) HC3N integrated intensity map. Plus symbols show the positions of the maser sites G305A and G305B. (b) HC3N spectrum 
towards G305A. The dotted line is a best fit model to the HC3N hyperfine structure, holding the optical depth constant at 0.1. Because 
we do not see any signs of the hyperfine components, we will assume the emission is optically thin. 
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Table 2. Properties of detections and non-detections. The notation for each Hne either conforms to the structure of (Juppcr - Jiowor) 
or (J_K,+K ^ J— k,+k)- Only detections for each transition are reported, except where no detection anywhere was made. In the case 
of a non-detection, then la upper limit for the integrated intensity is given. The last column lists names that are used in the text as a 
reference for the position of the emission. 



Observed Position Integrated Peak Line Association 

Line (J2000) Intensity Velocity Width 

(K.kms-i) (kms-i) (kms"!) 



HC3N (2-1) 




13 


11 


13 


71 


-62 


34 


40 


5 


CH3OH (32.1 -3l,2) 


E 


13 


11 


13 


71 


-62 


34 


41 





CH3OH (42,2 -4i,3) 


E 


13 


11 


13 


71 


-62 


34 


41 





CH3OH (22,0 - 2i,i) 


E 


13 


11 


13 


71 


-62 


34 


41 





CH3OH (32,1 -3l,2) 


E 


13 


11 


14 


00 


-62 


34 


46 


5 


CH3OH (42,2 -4i,3) 


E 


13 


11 


14 


00 


-62 


34 


46 


5 


CH3OH (22,0 - 2i,i) 


E 


13 


11 


14 


00 


-62 


34 


46 


5 


CH3OH (32,1 -3l,2) 


E 


















CH3OH (42,2 - 4i,3) 


E 


13 


11 


10 


89 


-62 


34 


39 





CH3OH (22,0 - 2i,i) 


E 


















NH3 (1,1) 




13 


11 


13 


71 


-62 


34 


41 





NH3 (2,2) 




13 


11 


13 


71 


-62 


34 


41 





OCS (2-1) 




13 


11 


14 


08 


-62 


34 


39 





NH3 (3,2) 




13 


11 


13 


71 


-62 


34 


41 





H2O (61,6 -52,3) 




13 


11 


13 


71 


-62 


34 


41 





H20 (61,6 - 52,3) 




13 


11 


13 


50 


-62 


34 


31 





H2O (61,6 - 52,3) 




13 


11 


08 


29 


-62 


34 


45 






CH3OCHO (2i,2-li,i) E 
CH3OCHO (2i,2-li,i) A 

CH3OH (82,6 - 81,7) E 13 11 10.89 -62 34 39.0 
DC3N (3-2) 
CH3CN (1-0) 
DC3N (2-1) 

H91/3 (RRL) 13 11 08.87 -62 34 41.0 

HC7N (15-14) 
HC5N (7-6) 



9.6 


-41.9 


5.1 


G305A 


12.2 


-39.8 


5.1 


G305A 


11.8 


-39.8 


4.8 


G305A 


9.0 


-40.0 


5.9 


G305A 


4.0 


-39.7 


2.5 


G305A(SE) 


5.3 


-39.7 


2.1 


G305A(SE) 


2.3 


-39.8 


2.9 


G305A(SE) 


<1.2 








46 


-42.1 


<0.75 


3" E of G305B 


<1.2 








84 






G305A 


11 






G305A 


3.9 


-42 


<12 


G305A 


15.0 


-39 


<13 


G305A 


18000 


-40 


<3.4 


G305A 


25 


-52 


<3.4 


10" N of G305A 


6.3 


+4 


<3.4 


14" SW of G305B 


<0.11 








<0.11 








130 


-43 


<12 


3" E of G305B 


<0.03 








<0.05 








<0.03 








3.7 


-43 


<18 


G305HII 


<0.03 








<0.11 









then these spectral types will be revised to 08 for G305HII 
and 09. 5 for G305HII(SE). Note that we use the value of 
90% absorption to illustrate how this affects the spectral 
type. We do not consider 90% as the most likely figure dust 
over, say 10% or 99%. 

We do not see any continuum emission associated with 
either of the methanol maser sites G3Q5A or G305B, with a 
3a detection limit of 0.15 mjy and 0.09 mjy, respectively. 
This is suprising as both sites are expected to harbour 
massive stars, which should have detectable radio contin- 
uum emission, as is the case for G305HII and G305HII(SE). 
G305A in particular appears to be the concentration for 
strong millimetre line emission (Walsh et al. 2006, and this 
work). This leads us to the conclusion that G305A at least, 
and probably G305B, are at a very early stage of evolu- 
tion where any Hii region has not had enough time to form 
and/or expand to an observable size and strength. 



4 DISCUSSION 

4.1 Maser emission in G305.2+0.2 

As well as the two Class II 6.7 GHz methanol maser sites 
G305A and G305B, a number of emission features detected 
in this work may also be masers. Water masers are detected 
at three positions in the field: coincident with G305A, ap- 



proximately 10" N of G305A and approximately 14" SW of 
G305B. This last maser site is coincident with the extended 
radio continuum emission of G305HII (Figure |3} . The three 
water maser features appear over a velocity range of about 
56 km s~^ (Figure [SJ. It is therefore likely that at least one of 
these features arises in an out flow, as water masers are com- 
monly found in outflows (eg. iMoscadelli. Cesaroni fc Rioial 
(2005)), which can give rise to such high velocities. 

The J=4 and J=8 methanol features exhibit very bright 
and narrow emission about 3" to the east of G305B, where 
the Class II 6.7 GHz methanol m aser lies. We suspect that 
both are Class I methanol masers (jMiiller. Menten fc Maderl 
[2004'). It is suprising that we see Class I and II methanol 
masers in close proximity to each other as it is unlikely that 
both Class I and II m asers can be strong in the same place 
jVoronkov et al.ll2005^ . However, we note that the 3" pro- 
jected offset between the two corresponds to 12000 AU (as- 
suming a distance of 3.9 kpc), which is certainly well sepa- 
rated enough that there may be different physical conditions 
at the positions of the Class I and II masers. Figure [1] shows 
that G305A(SE) is a narrow line emission feature in the J=4 
methanol line, which may also be a maser. Due to its intrin- 
sic weakness, we cannot rule out the possibility that it is not 
a maser, but unusually narrow-lined thermal emission. We 
explore the possibility that this is thermal emission in the 
next subsection. 
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(a) G305.S+0.g 



18496 MHz Continuum 



C 




o 



(b) G305.340.2 



Right Ascension (J2000) 



9 



G305HII(SE) 



Right Ascension (J2000) 

Figure 3. (a) 18496 MHz continuum emission towards G305. The 
two 6.67 GHz methanol maser sites G305A and G305B are indi- 
cated with the plus symbols. Water masers are indicated with 
diamond symbols. The contours are at 10, 20, ... 90% of the peak 
flux density, which is 22mjy/beam. (b) Spitzer 8fim image to- 
wards G305. Very strong emission is seen concident with G305B. 
This creates an artifact that extends to the NE of G305B. No 
emission is seen coincident with G305A. 



4.2 Thermal CH3OH (methanol) 

For G305A, integrated intensities of the three thermal 
methanol lines are listed in Table [2] Assuming the lines 
are optically thin, we can estimate the methanol column 
density a nd rotation temperature using a ro tation diagram 
analysis ijLinke. Frerking fc Thaddeuslll979h . shown in Fig- 
ure|6] We find a column density of (9.6 ±4) x lO^^cm"^ and 
temperature of 26 K. The errorbars shown in Figure [6] are 
determined by the relative uncertainty in the line integrated 
intensities. This uncertainty derives from the fit performed 
to each line, as well as the inherent uncertainty in the image 
reconstruction process. However, it does not include an abso- 
lute calibration error. We estimate that the error in our fiux 
determinations is approximately a factor of two, based on 
the flux calibration of the telescope. The size of this error is 
confirmed by comparing observations of the J=3 and 4 tran- 



C305HI1(SE) 



Frequency (GHz) 



Figure 4. Continuum flux as a function of frequency for G305H1I 
and G305H1I(SE). The general trend is a flat distribution, consis- 
tent with optically thin free-free emission. The scatter seen in the 
points is likely due to the errors in the absolute flux calibration 
for each frequency and is therefore a good representation of the 
error in the flux of each Hll region. 



40 
30 
20 




Velocity (km/s) 



Figure 5. Water maser spectra for the three detected features. 
The strongest feature is shown in the top spectrum and is coinci- 
dent with G305A. This feature also dominates the emission seen 
in the other spectra. The dotted lines in the middle and bottom 
spectra show a 10-fold increase in the intensity axis, to highlight 
the two weaker maser features that appear 10" to the north of 
G305A (middle) and 14" to the SW of G305B (bottom). 



ATCA Observations of G305. 2+0. 2 7 



42 



38 



42 



38 



4,, - 4 . 



Troi = 26K 

K = 9.6 xlQis cm- 



" 1 



N - LI xlO'^' cm-2 



G305A(SE) 



^26 ^17 



20 40 60 80 

Energy of upper state (K) 



100 



Figure 6. Methanol rotation diagram for G305A (top) and 
G305A(SE) (bottom). For each source, a straight (dashed) Une 
is fitted through the three lower temperature transitions (J=2,3 
and 4). The fitted rotational temperature and column density is 
given for each line. Whilst the J=8 upper limit is consistent with 
the fitted line for G305A, it is not possible to fit a line that is con- 
sistent with the J=8 upper limit for G305A(SE). The dotted line 
in the lower plot shows a fit through the J=8 line, which results in 
a colum 
of 78 K. 



4.3 HC3N (Cyanoacetylene) 

Cyanoacetylene possesses hyperfine structure in its (2- 
1) transition. However, no hyperfine components were 
observed, apart from the main two transitions at 
18196.279 MHz. Figure [2] shows a fit to the data, assum- 
ing a total optical depth of 0.1. The satellite components do 
not appear above the modelled spectrum. Therefore, we can 
be confident the main line cyanoacetylene emission has an 
optical depth less than 0.1 (ie. optically thin). 

Given the integrated intensity listed in Figure [21 and as- 
suming an excitation temperature of 26 K from the methanol 
measurements, we calculate the cyanoacetylene column den- 
sity is 4.0 X 10"cm"^. 



4.4 NH3 (Ammonia) 

Ammonia is detected towards G305A as an unresolved 
source. Because our low spectral resolution data cannot re- 
solve the hyperfine structure of the (1,1) and (2,2) inver- 
sion transitions, we are unable to use our data to inves- 
tigate the optical depth of the ammonia emission. How- 
ever, higher spectral and spatial resolution ammonia obser- 
vations (Steven Longmore, private communication 2007) de- 
tect (1,1), (2,2) and (4,4) emission coincident with G305A. 
The (1,1) and (2,2) transitions appear optically thick, so it 
is not possible to derive a reliable ammonia column density. 
The presence of (4,4) emission, which has an upper energy 
level of 200 K above ground, strongly suggests a hot compo- 
nent to G305A. 



sitions of methanol performed during both sessions in May 
and September. In addition to the three detected methanol 
lines, we also observed the 826 — 817 transition of methanol, 
but failed to detect it. Therefore, we include the upper limit 
for this transition in Figure [6] We find that the upper limit is 
consistent the best fit line to the other three points. Because 
the emission comes from regions which appear unresolved, 
we do not consider missing fiux due to missing short UV 
spacings a significant problem. 

For G305A(SE), we also fit the three methanol lines 
in Figure [6l as well as the non-detection in the J=8 line. 
Here we find a column density of 1.1 x lO^'^cm"^ and a 
temperature of 78 K. The derived excitation temperature 
from the rotation diagram analysis does not have a phys- 
ical interpretation as a kinetic temperature unless the gas 
density is above critical. Hence we perform, in i|4.8l a non- 
LTE analysis to derive column densities and relative abun- 
dances. As mentioned in M.ll the narrow linewidths of the 
detected methanol transitions at this position are sugges- 
tive of maser emission. In addition to this, the methanol 
spectrum at G305A(SE) is similar to the methanol maser 
spectrum 3" to the east of G305B - the J=4 feature is the 
strongest. Such maser emission could explain the unusual re- 
sults of the rotational diagram analysis for G305(SE) and so 
we favour the explanation that this is weak maser emission. 



4.5 OCS (Carbonyl sulfide) 

Carbonyl sulfide is detected as an unresolved source towards 
G305A. Assuming an excitation temperature of 26 K, we cal- 
culate a column density of 2.5 x lO^^cm"^ for this molecule. 



4.6 H91/3 radio recombination line 

This radio recombination line is detected in the field and 
appears coincident with the Hii region G305HII. 



4.7 G305SW - a possible prestellar massive core 

IWalsh fc BurtonI (|2006l) reported the detection of G305SW 
in N2H+, but found it lacking in detectable emission in other 
tracers. They suggested this may be the signature of an ex- 
tremely young and cold core, potentially prestellar. Unfor- 
tunately, G305SW lies just outside the primary beam of our 
2005 May observations, but lies within the primary beam 
for the 2005 September observations. We found no evidence 
for emission in any of the tracers from the 2005 Septem- 
ber observations. This is not suprising for most of the more 
complex molecules, but it is suprising that we did not de- 
tect any ammonia associated with G305SW. Since N2H"'" 
and ammonia are chemically very similar, we would expect 
them to trace similar regions. 
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4.8 Derivation of column densities 

The column densities listed above were derived with the as- 
sumption that the gas is in local thermodynamic equilib- 
rium (LTE). The rotation diagram method used (see for 
example Figure 5) produces inaccurate temperatures and 
column densities if only limited line data is available (<3 
lines). Moreover, the assumption of LTE is only valid for 
large densities. We present below, column densities derived 
from a non-LTE analysis of the molecular lines observed in 
Table 2. 

Observations have shown CH3OH accounts for a 
significant fraction of ice mantles relative to H2O 
iTielens fc Allamandolalll987l '). It is believed that CH3OH is 
formed on the surfac es of grains by the hyd rogenation of CO 
on the grain surface (|Charnlev et al.lll992f ). Above tempera- 
tures of ~ 90 — 100 K, all ice species, including CH3OH will 
be sublimated from the surfaces of grains. The abundance 
of CH3OH should therefore increase in the inner warmer re- 
gions of hot cores. In §4.2, three CII3OII lines were analyzed 
with a rotation diagram and a kinetic temperature of 26 K 
was derived. However, to release the methanol from the sur- 
faces of grains, the temperature must have been of the order 
90 — 100 K at some earlier time. It seems unlikely that the 
gas would have rapidly cooled to 26K, and we attribute the 
low 26 K temperature to the inaccurate assumption of LTE. 

In order to account for the discrepancy in the derived 
low (26K) temperature using the rotation diagram analysis 
of the CH3OH lines, we present below the results from a non- 
LTE radiative transfer model. A new analysis of the molecu- 
lar lines in Table 2 produces new column density values and 
temperatures, which agree better with our understanding 
of the underlying physics and chem istry (see also 35|. For 
this calculation, molecular data from lSchoier et ahl (|2005h is 
used which includes energy levels, statistical weights, tran- 
sition frequencies as well as coUisional rate coefficients (for 
collisions with II2), for each spec ies of interest. In partic- 
ular, we use CH3OH rat es from iPottage. Flower fc Davia 
||2004). NH.3 rates from iDanbvetaLl lll988h. OC S rates 
from breen fc ChapmanI (|l978l) and ISchoier et al] (|2005l ) 
and HC3N rates from Green fc Chapm"anl ( 1978h . The cal- 
culation is similar to that of Schoier et al.l (|2005l ). and cal- 



Integrated Intensity (K km/s) 



culates the strength of the emission line for a species and 
transition of interest, for a given density, temperature and 
column density. 

To derive a column density for a given species, the inte- 
grated fine intensity, J T^dv in K.km/s, is found as a func- 
tion of both density (uhj ) and column density. The source 
is assumed to fill the beam. Initially we assume some tem- 
perature and the line width for each line of interest is taken 
from the observations in Table 2. The constant contour of 
the observed integrated line intensity can then be traced. 
This will usually give an upper or lower limit to the column 
density. With more than one line, it is often possible to si- 
multaneously derive both a column density and density, for 
a given temperature. 

Figure [7^ shows a plot of the integrated intensity as a 
function of and column density assuming a temperature 
of 26 K for the CH30H(42,2 - 4i,3)E line from Table 2, with 
a FWHM of 4.8kms"\ The observed 11.8 K km s"^ con- 
tour is plotted as a dashed line. As the density is increased 
and LTE approached, the intensity becomes constant with 
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Figure 7. Contours of integrated intensity for rotational 
methanol transitions, (a) Integrated intensity of tlie CH3OH 42,2— 
4i 3 transition calculated using a non-LTE theoretical model as- 
suming a temperature of T=26 K. The constant contour which 
describes the observed intensity from Table [2] is plotted as the 
dashed line, (b) Comparison of the CH3OH 42,2-4i,3 observed 
intensity from panel (a) with the CH3OH 22,o^2i 1 and 32,i-3i,2 
transitions calculated in a similar way for T=26 K using the val- 
ues from Table [2] (c) and (d) are similar to panel (b), but for 
temperatures of 100 K and 200 K, respectively. 
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Table 3. Observations G305A from Steven Longmore (private 
communication 2007). 



Observed line 


Integrated Intensity 


Line Width 




(Kkms-i) 


(kms"-*-) 


NH3(1,1) 


33.58 


5.5 


NH3(2,2) 


23.05 


5.1 


NH3(4,4) 


22.23 


7.9 



column density. iHill et al.l ()2005l l report a density of the or- 
der ~ 1 — 2 X lO'^cm"^, implying a column density of 
~ 4 X lO^^cm"^ from Figure [7^. However, the two other 
lines observed for CH3OH can be used to further restrict 
the column density. 

In Figure [TJa, the corresponding constant contours of 
each observed integrated line intensity for all three CH3OH 
lines are over-plotted. All three lines overlap and agree for 
> lO^cm"'^, yielding the L TE limit of the column 
density. Since the observations of iHill et al.l l|2005h suggest 
a lower density, it is likely that the gas is not in LTE. 
In Figures [7}; and [TJi, the same lines are over-plotted as- 
suming temperatures of lOOK and 200K. In each of these 
cases, there is agreeement for densities and column densi- 
ties in the ranges 1.8 x 10*cm~'^ < < 8 x 10*cm~^ and 
7.9 x W^^cm-^< N <2.0x 10^^ cm"^ respectively. Since all 
three lines agree for both temperatures, this method does 
not constrain the temperature. However, as previously men- 
tioned, temperatures of 90 K or more are expected to evap- 
orate CH3OH off grain surfaces, and so we favour a temper- 
ature of at least 100 K. 

High spectral and spatial resolution ammonia observa- 
tions have been carried out (Steven Longmore, private com- 
munication 2007) for the (1,1), (2,2) and (4,4) lines. The 
details of the detection of these lines coincident with G305A 
are given in Table [S] A similar analysis to that of Figure [7] 
for the three NH3 lines is shown in Figure |8l In Figure [8^, 
the temperature is 26 K, and clearly there is no agreement 
between the NH3(4,4) line and the other two lines. However, 
if the temperature is increased, as in Figures [Sja and|8j; then 
there is closer agreement. The value derived from these plots 
for the column density is a little ambiguous since the (1,1) 
and (2,2) contours only converge as decreases, so we 
only quote a lower limit of Ntqug > lO^^cm"^. However, it 
is important to note that the observation of the NH3(4,4) 
strongly suggests that there is a temperature in excess of 
100 K. 

We note that densities of 1.8 x 10* - 8 x 10*cm"^ 
are unusually low for traditional hot core densities (~ 
lO^'cm"'^). However, complex organic molecules are seen 
in extended emis sion in Sgr B2 and in the central molecu- 
lar zo ne (Reguen a- Torres et al. 2006). Requcna- Torres ct al. 
I2OO6) find that these molecules, including CH3OH, can have 
low excitation temperatures (10-20 K) despite high kinetic 
temperatures (~100K) and cold dust (<20K) because the 
density is low (~ 10*cm~'^). They argue that shocks remove 
mantle material and drive 'hot core' type chemistry. So low 
excitation temperatures may not necessarily reflect cold gas- 
phase chemistry. 

We cannot do a complete analysis of OCS and HC3N, 
with only single transitions for each. However we can give 
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Figure 8. Contours of integrated intensity for rotational am- 
monia transitions, (a) Integrated intensities of tlie NH3 rotation 
transitions (1,1), (2,2) and (4,4) calculated using a non-LTE the- 
oretical model assuming a temperature of T=26K. The chosen 
contours are taken from the observed intensities of Table \3\ In 
panels (b) and (c) the same comparisons are made, assuming 
temperatures of 100 K and 200 K, respectively. 



new values assuming densities in the range 1.8 x lO^'cm"^ < 
< 8 X 10*cm"^. Figure H shows the OCS(2-l) and 
HC3N(2-1) lines plotted for lOOK and 200K. For densities 
in the range 1.8 x lO'^cm"'^ < riHa < 8 x lO'^cm"'', we can 
put limits on the column densities, which are given in Table 
ll for lOOK, along with the other non-LTE values for NH3 
and CH3OH. We note that column densities derived here 
compare well to those derived in sections 14.21 14.31 and 14.51 
which assume LTE. This indicates that column densities de- 
rived under the assumption of LTE, at least in the case of 
G305A, may provide reasonable estimates of the true col- 
umn density. 
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Figure 9. Contours of integrated intensity for rotational OCS 
and HC3N transitions, (a) Integrated intensities of tlie OCS (2- 
1) and HC3N rotation transitions calculated using a non-LTE 
theoretical model assuming a temperature of T=100K. The cho- 
sen contours are taken from the observed intensities of Table [2] 
(b) the same transitions are presented as in panel (a), but with a 
temperature of 200 K. 



Table 4. Column densities for G305A from igS] for 1.8 X 



lO^cm-^ < nn, < 8 x lO^cm" 



Molecule 



Column Density (cm ^) 



CH3OH 
NH* 

OCS(IOOK) 
OCS{200K) 
HC3N(100K) 
HC3N(200K) 



'<N < 2.0 X lO^^'cm- 



7.9 X lO^^c 
> 10^5 

4.0 X 10^5 < ^Q^g < g 3 X iQis 
5.6 X 10^5 < ATqcs < 9-5 X 10^5 
2.0 X 10" < AfHCaN < 6.3 X lO^"' 
2.8 X 10^4 < AfHCaN < 6.3 X lOi* 



'value derived from observations of Steven Longmore (private 
communication, 2007). 



5 TIME-DEPENDENT CHEMISTRY 

The analysis in section 4, in particular the detection of the 
NH3 (4,4) line, strongly suggest that G305A has a hot com- 
ponent with a temperature of the order 100 - 200 K. We 
will therefore now model the chemistry of G305A, based on 
the assumption that it is a hot core. The time dependent 
chemistry of a spherical core, representative of G305A, will 
be calculated below and the observations will be discussed 
in hght of these results. A molecule found inside a hot core 
may be present due the following three processes: if it is ei- 
ther formed by cold chemistry in the earlier core collapse 
phase then accreted onto a grain and re-released inside the 
hot core, formed by grain surface chemistry and released, or 
formed in the gas phase chemistry of the hot core. Although 
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Figure 10. Time dependent column densities for the four species 
of interest, as calculated from a spherical hot core chemical model 
of G305A as described in [[5] A constant density is assumed, 
njjj = 1-8 X 10'*cm~'^, with a temperature of 100 K (top panel) or 
200 K ( bottom panel). A comparison of a species column density 
constrains the age of the core. The grey shaded bar illustrates the 
age range of the core derived in comparing the observed column 
densities given in Table |4] 



the mantle compositions are relatively simple e.g. CH3OH 
and NH3 , the release of these into the gas phase are believed 
to be the drivers of the very complex chemistry observed in 
hot cores. The relative abundances of the mantle species, 
along with varying densities and temperatures can lead to 
chemical differences among cores. 

We now attempt to model the gas phase chemistry of 
G305A. We assume the core is spherical with a density 

= 1-8 X lO'^cm"'' and uyi^ = 8 x 10^cm~^ based on 
the a nalysis in section f48l and a radius of 1.1 pc (|Hill et al.l 
l2005h . The model is run twice with constant temperatures 
throughout o f 100 K and 200 K. We use a model similar to 
iMillar et al.l (|l99'i1 ). which consists of 22 depth points, or 
shells, each with a specified density and temperature. Col- 
umn densities are calculated by integrating along lines of 
sight, and to account for beam dilution, the column densi- 
ties are weighted and summed oyer a gaussian beam (eg., 
[Thompson. Macdonald fc Millad Il999l ). The source is as- 
sumed to fill the beam. 211 species are chosen with 2575 
reactions and rate coefficients, taken from the online UMIST 
database l|Le Teuff et al.ll2000l ). 

The ini t ial ab undances are taken from Table 2 of 
[Miliar et al.l l|l997l ) which assumes that the gas is mostly 
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Figure 11. Column densities from tlie splierical core model, 
similar to Figure 1101 but with the density increased to riHa = 
8 X lO^cm-^. 



molecular with some ionization by cosmic rays. Assum- 
ing spontaneous evaporation from grain surfaces, there are 
also the primary injected mantle molecules CH3OH, NH3, 
C2H2,C2H4,C2H6,CH4, H2CO, CO2, O2, H2O, H2S. Ice ab- 
sorbtion spectra towards massive young stellar objects sug- 
gest that OCS sho uld also be included as a mantle species 
iKeane et aTll2001^ . We also find that an initial abundance 
of OCS is needed in the model to generate an OCS abun- 
dance that resembles the observations, the gas phase chem- 
istry alone cannot produce the observed column density (see 
the discussion for Figures[lO]and[Tl]below) . The initial abun- 
dances are g iven in Table[5] which are based o n the models of 
iMillar et all ([1997) and No mura and Millar! (|2004 ). The ini- 
tial conditions are also chosen so that the total abundance of 
C, N and O do not exceed with the mean interstellar values; 
C/H = 1.4 X IQ-^, 0/H= 3.19 x 10"^ and N/H = 7.5 x 10"^ 
l|CardeUi et al.lll996l : [Meyer et al.lll997l . Il998l ). It should be 
noted that differences in the mantle chemistry and choice of 
initial abundances could drive a different gas phase chem- 
istry, this will be considered in future papers, and is outside 
the scope of the current work. 

The model results are presented in Figures [1111 and lTTI for 
riHa = 1-8 X 10*cm~'^ and = 8 x lO^cm"^, respectively. 
The column densities of the four species of interest are plot- 
ted as a function of time since the onset of the core heating. 
In each figure the chemical model has been run twice as- 
suming temperatures of 100 K and 200 K. By comparing the 



Table 5. The initial abundances, with respect to H nuclei, assumed 
in the chemical model 



Species 


Abundance 


Species 


Abundance 


CH3OH 


5 x 10"'^ 


C2H2 


5 X 10"^ 


NH3 


1 x lO-'^ 


C2H4 


5 X 10"'^ 


H2S 


1 X 10-6 


C2H6 


5 X 10-9 


H2O 


1 X 10-5 


CH4 


2 X lO-'^ 


O2 


1 X 10-6 


CO2 


5.0 X 10-6 


H2CO 


4 X 10-* 


OCS 


5.0 X 10-* 


CO 


5 X 


N2 


2 X 10-'' 


He+ 


2.5 X IQ-" 


H+ 


1 X 10-1° 




1 X 10"* 


Si+ 


3.6 X 10-* 


Fc+ 


2.4 X 10-* 







column densities for each species in Tabled derived from the 
observed emission lines in Tabled the chemical model can 
be used to constrain both the temperature and age of the 
core. In th e gas phase, the produ ction of CH3OH is very inef- 
ficient (eg. lWoodall et all j20o3)); g oing via radiative associ- 
ation of CHj " and H2O at a l ow rat e ([Luca. Voulet fc GerlichI 
I2002D . while [Ceppart et al.l (|2006l ) have recently shown that 
the fraction of dissociative recombinations which produce 
methanol are of the order 3%. From the methanol column 
density limits quoted in Table O upper age limits are de- 
rived in Figure [TTI for 71^2 = 8 x 10*cm-^ ; 3.1 x 10"* years 
for 100 K and 2 x 10"* years from 200 K. The lower density 

= 1-8 X lO^cm-'^ cases of Figure [Tlll do little to constrain 
the age since Neuron < 7.9 x lO^^cm-^. From the NH3 ob- 
servations, both Figures [Till and [TTI give upper age limits of 
1.6 X 10^ years. 

The column density of HC3N indicates either, a mini- 
mum age of r-^ 10'* years, or an age of 10^ years. We favor the 
later time b ased on the NH3 co mparisons discussed above. 
As noted bv lMillar et al.1 l| 19971 ). HC3N is usually regarded 
as a cold cloud tracer, however it is formed in a hot gas 
through reactions of N and C3H2 and C2H2 with CN, where 
C2H2 is evaporated from grain mantles. This suggests that 
HC3N should be a useful species to trace the conditions and 
ages of hot cores. The highly variable HC3N abundance with 
time, in Figures [TUI and [TTI further supports the use of HC3N 
as a 'chemical clock' species. 

The OCS column densities quoted in Table [4] are > 
4 X 10*5 cm-"^, which is higher than the maximum OCS col- 
umn density for all cases. The observed OCS abundance can 
be reproduced in the model if the initial abundance of OCS 
was increa sed. For simplicity, we ass umed the OCS abun- 
dance from lNomura and Millaij (|2004l ) based on their model 
of G34. 3-1-0. 15, and so it is highly possible that the initial 
OCS abundance was greater for G305A. We will examine the 
effects of changes in the model for varying initial conditions 
in future papers. It is important to note that if no initial 
abundance of OCS is assumed in the model, then the gas 
phase chemistry only produces a OCS column density of the 
order 4 x lO^^cm-'^ (plots not shown here), which strongly 
suggests that the origin of the OCS is the evaporated ice 
mantles. Moreover, OCS has be en detected in ice- absorbtion 
features of massive protostars ([Keane et al.ll200l[ ). 

It is important to note that the column densities of 
NII3 and IIC3N give similar core ages, lO^years, for each 
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of the 4 cases plotted in Figures [10] and [TT] However, the 
determination of the temperature and density is still a little 
ambiguous. There is better agreement for NH3 and HC3N 
times for the denser hotter case. Moreover, the NH3 ob- 
servations (Steven Longmore, private communication 2007) 
discussed in H4.8I ( in particular Figure [S} suggest higher 
temperatures of the order 200 K. We therefore favor the high 
density 8 x lO'^cm"'' and 200 K case, suggesting a core age 
in the range 2 x 10"* < tcore < 1-5 x 10^ years. 
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6 CONCLUSIONS 

We have observed the massive star forming region 
G305.24-0.2 at 1.2 cm. We detected emission in methanol, 
cyanoacetylene, ammonia, carbonyl sulfide and water, as 
well as the radio recombination line H91/3. All molecular 
transitions are detected towards G305A, which confirms the 
hot core nature of this source. Three lines of methanol are 
found in emission at the position of G305A, and together 
with a fourth transition which was not detected, we have 
constructed a rotational diagram where we derive a rota- 
tional temperature of 26 K and a column density of 9.6 
X 10^^ cm~'^ for methanol. However, because methanol re- 
quires temperatures in excess of 90 K to evaporate from dust 
grains, we find the temperature unrealistic and attribute this 
to the incorrect assumption of LTE. With a non-LTE anal- 
ysis of G305A, we derive a model which is consistent with 
a density of 8 x lO'^cm"^, and a temperature of the order 
200 K. A gas phase time dependent chemical model further 
constrains the age of the core to 2 x 10* < tcor& < 1.5 x 10^ 
years. 

Extended emission to the SE of G305A is seen in three 
methanol transitions. It has a narrow line width and an 
unusually high rotational temperature (78 K). We interpret 
this as evidence that the emission may be a weak maser. 

Strong methanol maser emission is found 3" to the east 
of G305B. Since the methanol transition at 24.933 GHz is a 
Class I maser, G305B appears to be a unique object wor- 
thy of further investigation as it shows both Class I and II 
methanol maser emission within a small region. No other 
detections were made at the position of G305B. Since it is 
also known to be a luminous source, the lack of molecular 
emission suggests G305B is an older source that has already 
had time to clear out its surrounding molecular material. 

Radio continuum emission at 18.496 GHz is detected 
towards two Hii regions (G305HII and G305HII(SE)). 
G305HII shows a cometary morphology and is proba- 
bly powered by a star in the range 08 to BO, whereas 
G305HII(SE) shows a shell morphology and is probably 
powered by a star in the range 09.5 to BO. 

No radio continuum emission was detected towards ei- 
ther G305A or G305B. For the case of G305A, this may be 
because the source is too young to have produced an observ- 
able Hii region. However, it is not clear why the older, and 
infrared-luminous, source G305B does not show any radio 
continuum emission. 
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